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Abstract 
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The productions of muon pairs from the decay of heavy quarkonia have 

O . been evaluated for different centrality of the nuclear collisions at LHC energies. 

| The effects of the various comover scenarios on the survival probability of the 

' heavy quarkonia have been considered. The effects of shadowing and comover 

. suppressions on the dilepton spectra originating from the decays of J/tp is 

\ found to be substantial. The dilepton yield from the thermal J/ip has also 

J> ' been estimated and found to be small. 

: 

o ■ 

I. Introduction 

Ever since the possibility of creating quark gluon plasma (QGP) in relativistic heavy 
ion collision was envisaged, numerous signals were proposed to probe the properties 
of such an exotic state of matter. In this context Satz and Matsui [1] had suggested 
that the production of heavy quark resonances (J/ip) will be suppressed as a result of 



colour Debye screening in a hot and dense system of quarks, anti-quarks and gluons. 
This suppression could be detected experimentally through the dileptonic decay 
mode of these resonances. ALICE dimuon spectrometer [2] is dedicated to look for 
this type of signal. However, it is a daunting task to disentangle the contributions of 
the heavy quarkonium states to muon spectrum due to the background from several 
other sources, e.g. Drell-Yan, semileptonic decay of open heavy flavoured mesons 
(DD, BB) etc. Low energy muons from kaons and pions also constitute a large 
background. 

In this work we shall estimate the dimuon production from the decay of both 
hard {i.e. J/ip produced from initial hard process, will be called hard J/ip hereafter) 
and thermal J/i/)'s. In heavy ion collisions the production and decay of hard J/ip 
or T proceeds through the following three steps: (i) the production of pair of heavy 
1 Part of the work has been done at CERN Theory Division, CERN, Geneva 
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quarks (perturbative), (ii) their resonance interactions to form the bound state (non- 
perturbative) and (iii) the propagation of the quarkonia through the medium and 
their subsequent decay to dileptonic modes with certain branching ratios. 

The initial state in relativistic heavy ion collisions consists of either hadronic 
matter or QGP depending on the incident energies of the colliding nuclei. At LHC 
energies the formation of QGP is almost unavoidable. Even if the system is formed in 
QGP phase it will revert back to hadronic phase due to the cooling of the expanding 
QGP system and hence the interaction of the J/ip formed in initial hard collision 
with the hadronic matter is inevitable. Therefore,, we need to consider the survival 
probability of those J/ip due to its interactions with the hadronic medium. 

At high temperatures a non-negligible number of J / ip mesons is expected in the 
thermal medium, which decay to lepton pair. The dilepton spectra originating from 
these thermal J/ip is also estimated here. 

The paper is organized as follows. In section II we shall describe the formalism 
for the production, propagation and decay of hard and thermal J/^'s. In section 
III results of our calculations will be presented followed by summary and discussion 
in section IV. 

II. Production, Propagation and Decay of Hard 
and Thermal J/ip 

a. Production 

In this section we shall consider the hard J/ip production in the colour evaporation 
model (CEM) [3] and their decays to lepton pairs. As mentioned before this consists 
of two stages: (i) production of a cc pair (perturbative process) and (ii) subsequent 
non-perturbative evolution into asymptotic states. We have considered those hard 
processes which can contribute to cc productions irrespective of their colour and 
spin-parity. The colour neutralization occurs by the interactions (one or more soft 
gluon emission) with the surrounding colour fields and this step is considered to be 
non-perturbative. In CEM quarkonium production is treated identically to open 
heavy flavour production with the exception that in the case of quarkonium, the 
invariant mass of the heavy quark pair is restricted below the open meson threshold, 
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which is twice the mass of the lowest meson mass that can be formed with the heavy 
quark. Depending on the quantum numbers of the initial QQ pair and the final state 
quarkonium, a different matrix element is needed for the resonance production. The 
effects of these non-perturbative matrix elements are combined into the universal 
factor F[nJ PC ] which is a process and kinematics independent quantity [4]. It 
describes the probability that the QQ pair forms a quarkonium of given spin (J), 
parity (P) and charge conjugation (C). The production cross section for a J/ip or 
T is therefore given by [4] 

a[R(nJ PC )]=F[nJ PC ]a[QQ], (1) 

where the non-perturbative (long distance) factor can be written in terms of the 
probability to have colour singlet state (1/9) and the fraction p R of each specific 
charmonium state. The perturbative contribution (short distance) is given by 

a[QQ]= L Q dM ^itM^- (2) 

The contributions to heavy quark production in leading order come from qq — > 
QQ and gg — > QQ. The differential cross-section for heavy flavour production in 
hadron-hadron collision is given by [5] 

da r , , nnvi H(x a ,x b ,Q 2 ) 

[h A h B -> QQX\ = , (3) 



dM 2 Q dy 



where 



H{x ai x b) Q 2 ) = Qf {xa,Q 2 W (x b .Q 2 )+ q h f A (x a ,Q 2 )q h f B (x b: Q 2 )a q ^ QQ 

f 1 J 
+g hA (x a , Q 2 )g hB (x bl Q 2 )a gg ^ QQ . (4) 

x a>b = Me^/^/s, ^/i being the centre of mass energy of the hadronic system, q/s 
and g's are the parton distribution functions (PDF) to be taken from CTEQ, MRST 
or GRV [6]. Combining eqs.(l),(2) and (3) we obtain the cross-section for resonance 
production per unit rapidity as, 

We take p R = 0.5 (0.207) for J /if; (T). 



Next we consider J/ip(T) production in p — A and A — B collisions. To this end 
we first briefly mention the necessary formulae in Glauber model [7, 8] The total 
inelastic cross-section in A — B collisions at an impact parameter b is given by 



1 - 



l-T AB {b)a{ 



NN 



-i AB 



= 1-P (b), (6) 



db 

where Tab (b) is the nuclear overlap function given by 

T AB (b) = J d 2 sT A (b)T B (\b-s\) (7) 

The nuclear thickness functions are normalized to unity, i. e. J d 2 bT A (b) = 
f d 2 bdz p A (b, z) = 1 

Generally we are interested in the cross-sections for a set of events in a given 
centrality range defined by the trigger settings. Centrality selection corresponds to 
a cut on the impact parameter b of the collisions. The sample of events in a given 
centrality range < b < b m , contains a fraction of the total inelastic cross-section. 
This fraction is defined by 

Jo 6 ™ db da> " 



/(U = 1 „-£b ( 8 ) 



Jfdb 



Now we discuss the J/ip survival probability when it propagates through the 
hot hadronic medium produced in nucleus- nucleus collisions. After creation the J/ ip 
meson can interact with other nucleons in the target and the projectile and may 
get destroyed mainly due to J ftp — N interactions. The cross-section for J/ijj(T) 
production in p — A collisions can be written as 

rbr, 



a j/4,( b m) = A dbdz p(b,z) exp -(A - 1) J a ahs p(b, z) dz' 



where a™ is obtained from eq.(5). The interpretation of the above equation is 
as follows. The resonance is formed at r — (b, z) where the density of the target 
nucleus is p(r). It can travel in forward direction (z) at constant impact parameter 
and its intensity is attenuated due to J/ip — N inelastic collisions. The exponential 
factor accounts for this attenuation loss. 

The generalization of eq. (9) in nucleus- nucleus collisions is straightforward. The 
J/ip(T) production cross-section in A — B collisions in the impact parameter range 



< b < b m can be written as 



drr AB drr NN 



'%/// ^ = — J <h) T) AB I dsd^d^PAis^Zx) p B {\b - s\,z 2 ) 
x exp (^~(A - 1) J o- abs p A (s, z) dz'^j 
x exp (-(5-1) f° a^ s p B (\b-s\,z')dz') (10) 



Z'2 



where <J^/ T n is obtained using the nuclear parton distribution functions (PDF) 
containing shadowing effects. 

b. Propagation 

The experimental data suggest that besides the absorption due to J/^-nucleon 
interactions, there are additional sources of J/ip absorption in nuclear collisions. 
The produced J/ip particles can interact with the relatively heavier mesons (such as 
p, i] etc.) expected to be produced at a proper time r . Such interactions can lead to 
the disappearance of J/tp. Not all the hadrons produced in nucleus-nucleus collisions 
(co)move with the J/ip. So the number density of comoving hadrons is given by 
^co = fc^h, where is density of produced hadrons and f c is the fraction that 
(co)moves with the J/ip . In the following sections we shall discuss three different 
'comoving' scenarios. 

Scenario I 

In this section we outline the comoving scenario of Ref. [8]. We consider that a 
hadronic matter is formed at a proper time r with density n/ l (r ) and evolves 
hydrodynamically. If we assume Bjorken's scaling [10] solution, then at a later time 
(r) the density of hadron (rihij)) is given by 

n h {T) = (11) 

During the time span from r to r F (the freeze-out time) the J / ip particles interact 
with the comoving fluid. The survival probability S co is given by [11] 

S co = exp[- dT(a co v rcl f c n h (T)}] (12) 

J TO 
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Using eqs.(ll) and (12) we get 



S co = exp[-( a co f re i ) r f c n h {r ) In 



?1 

TO 



(13) 



where <7 co is the J/^-comover cross-section leading to the breakup of J/ip particle, 
v Te \ is the relative velocity between the J/ip and the comover. We have used a 
constant value for a co in our calculation. 

The hadron density rih can be related to the multiplicity of produced hadrons 
at r from the row-on-row collisions of the projectile and the target nucleons. The 
volume of the row is a™ Tody at proper time r . The hadron density of this row is 
therefore: 



n h (T ,b,s) 



dN 



AB 



0in r o dy 



(14) 



where dNAB/dy(b, s) is the hadron rapidity density which can be calculated from 
the participant density, n' w (b, s). For a row-on-row collision, it is given by 



nUM=<C AT A (b)+BT B (b-s) 



and the multiplicity is obtained as 

dNAB t r -k 1 dNiyN NN 
-(0, Sj ~ - ; (7;, 



2 rfy 



AT A (6) + 5T B (6-s) 



(15) 



(16) 



where is the multiplicity in nucleon-nucleon collisions. Combining eqs.(12)- 

(16) we obtain the J ftp survival probability, S co in a comoving scenario as 

S co (b,s) =exp(-c [AT A (b)+BT B (b-s)]) 7 (17) 

where Co is given by 



r F ldN NN 
c = (o"co ^rci ) In — / c 



(18) 



t 2 

For hadron rapidity density in nucleon-nucleon collision we use the following form [13]: 
dN NN /dy = 2.5 - 0.25 ln(s) + 0.023 ln 2 (s). 

Scenario II 

In Ref. [5] it is assumed that all the comoving hadrons are generated from the 
participating nucleons, which might be a reasonable assumption at SPS energies. 
However, at higher energies it needs modifications, as discussed below. One would 
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expect that the total multiplicity in nucleus-nucleus collisions at RHIC and LHC 
energies will come from both soft (~ N paxt ) as well as hard (~ iV co ii) collisions. It 
has been shown recently in [13] that about 10 % of the total multiplicity at RHIC 
energies comes from hard collisions, i.e. the total hadron multiplicity in nucleus- 
nucleus collisions at a given b can be written as 



dN AB 
dy 



(b) 



;i _ x) + x Ncoii(b) 



NN 



dy 



(19) 



where 



N pSLr t(b) = J dsn w (b,s) 
N coll (b) = AB<j™T AB (b) 

The participant density n w (b, s) in this case is given by 

n w (b, s) = ATa(s) (l - [l - cr^T B (b- § 



(20) 



B 



+ BT B (b - s) {1 - [1 - a* N T A (s) 
The initial density of the comoving hadrons in this scenario is: 

^(ro)=n h (r )f c = f/^-^-. 

Thus the comoving survival probability is given by 

c fu\ \ i \ dN AB {b) f c t f 

S co {b) = exp -(a co w rel ) In — , 

dy nR A t 



(21) 



(22) 



(23) 



Scenario III 

We assume that the comoving density is proportional to the final dN A B / dy which is a 
function of centrality. In this scenario it is also assumed that the final multiplicity of 
the produced hadrons depends on both {N part ) and (N co \\) (two component model). 
In this scenario S co can be estimated as [5] 



S co (b, s) = exp[-( o- co v rc i ) n w (b, s) In — ], 

to 



(24) 



where it has been assumed that the number of comoving hadrons is proportional 
to the number of participants. At low energies (e.g. SPS) this might be a valid 



assumption. However, at higher energies a certain fraction of comoving hadrons 
may come from the hard collisions (oc N co \\). To take into account this fact we 
modify the above equation to obtain 



S<x>(b, s) = exp < - ( a co v rei ] 



n w (b, s) -> 
1 — x) V xn c (b, s) 



In - }, (25) 



TO 

where n w (b, s) is given by eq.(21) and n c (b, s) = AB af^ T A {b) T B {b — s). 

In relativistic heavy ion collisions the yield is the relevant quantity rather than 
the cross-section. In order to obtain the differential number distribution from the 
cross-section one has to resort to Glauber model of nucleus-nucleus scattering [8]. 
Incorporating the comoving survival probability in eq.(10) we obtain the J/ip pro- 
duction cross-section in A — B collisions as 

da AB da NN 

x exp ^- (A - 1) o- abs p A (s, z') dz'^j 

x exp (-(B - I) <j abs p B (\b- s\,z')dz' 
da NN 

= ABT{b m ) J j V ' (T) (26) 
dy 

The total number of J/ip(T) produced in A — B collisions can now be written 

as 

dy {bm) ~ ai B (b m ) dy [27) 
c. Decay of J/ip to lepton pairs 

To calculate h A h B — > RX — > l + l~ X from hard processes we proceed as follows. 
Symbolically we can write 

^[h A h B -> R(l + r)X] = V -^- d ^[h A h B -> RX] (28) 
dy T R dy 

After production the quarkonia propagate in the medium before decaying to 
lepton pairs with certain branching ratio. The finite width of the quarkonia may 
be taken into account by folding Eq. 28 with the spectral function (Ar(M)) of the 
quarkonia, 

^-[h A h B - R{l + r)X] = I dM 2 ^[h A h B - RX]^^-A R (M) (29) 
dy J dy I p 



where A R (M) is given by, 

1 MTr 

Using the above set of equations we obtain the differential cross-section for the 
lepton pair production from hard heavy quark resonance decay as 

-\h A h L 

7i yivi- - nr R r -i- ivi-L R 

n „ r2m nlR _ H(x„.Xh. M 2 ^) 

(31) 



1 






TT 


(M 2 


-m 2 R ) 2 




PR 


r2m D / B 


X 




1 I 




IJ 


2mn 



H(x a ,x b ,M 2 Q ) 
s 

Now the dilepton yield from hard J/ip(T) decay in a given centrality range 
< b < b m can be written as 

dN r ,„ 1 M 1 ^,- </5( T) 



- ^ (/+r)x] ^ (m 2 - M^fi d ; { {bm) (32) 

In a similar way one can also calculate the the total number of heavy quark 
pairs, in a given centrality range. These are: 



N^(b m ) = no™ {r) (33) 



where TZ is given by 

K 



(-^binary) ipm) 



a NN 

^in 



f»™ dbABT AB (b) 



^ db[l-{l-T AB (b)< N } AB ] 
d. Thermal J/ip Decay 



(34) 



The muon yields from thermal J ftp decay during the life time of the fire ball is 
expected to be small due two reasons: (i) less abundance of the J / ip in the thermal 
system of temperature ~ of few hundred MeV and (ii) the probability of the J / ip to 
decay within the lifetime of the hadronic system is a very small due to tiny width of 
the J/ip. However, for the sake of completeness we give the thermal spectra below: 
dN 2 J + 1 



-KR 2 A ZTf^MVj,^ l+l - 



dM 2 dy (2tt) 2 

Mr tot 



i 

x — 

TT 



(M 2 - m 2 ) 2 + (MT tot ) 2 



T i dT 



J/i>> 

roo rr]2 \ 

x zdz — (35) 

Jm/t J Vl exp{z cosh(y — rj)) — 1 
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Figure 1: The centrality dependence of the suppression factor, 
where z = Mt/T. 

As mentioned above due to the small width of the J / ip most of them will decay 
after the fireball freeze-out. The yield from these J/ip have also been estimated using 
Cooper-Frye formula [14] and added to the thermal contributions. The expression 
for this contribution is given below: 

dN _ gR 2 A r F 1 
dM 2 dy ~ 4tt 2 n 

oo 

x E(+) n+1 / K l (nm T /T)m T d 2 p T (36) 

n=l J 

III. Results 

The variation of J-(b m ), which is the ratio of J/ip production cross-section in A — B 
collisions (scaled to p — p) to that in nucleon-nucleon collision as a function of 
centrality (b m ) is shown in fig. (1). Here we have assumed that there is no suppression 
due to comoving hadrons. In presence of comovers J r (b m ) will decrease further. 
However, in the evaluation of dilepton yield we have incorporated the suppression 
due to comover in three different scenarios as discussed earlier. 

In fig. (2) we show the centrality dependence of the total J ftp production cross- 
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Figure 2: Total J/ip cross-section in different comoving scenarios (see text) as func- 
tion of centrality at LHC energies. 

section in nucleus-nucleus collisions. It is seen that the survival probability is very 
sensitive to the choice of ctj/^n cross-section as well as to the comoving scenario 
adopted. 

In the second scenario of comoving absorption the survival probability becomes 
function of b only. So in fig. (3) we plot S co as a function of impact parameter. 
The suppression depends on two factors: (i) what fraction (x) of produced hadrons 
come from hard processes and (ii) what fraction (/ c ) of these hadrons (co)move 
with the J/ip. To show the sensitivity on these two factors we have chosen various 
combinations of x and f c . It is seen that for a given f c as x increases (which is 
possible as the beam energy increases) the survival probability goes down. 

Now let us turn to the lepton pair yield from J / ip decay. In fig. (4) we compare 
the lepton pair yield from Pb — Pb collisions at y/s = 6000 GeV for - 10 % 
centrality. We have chosen different combinations of x and f c to show the sensitivity 
of the results in second scenario of comoving suppression. However, we have checked 
that for very small values of x and f c the three scenarios for comoving suppression 
coincide. But for higher value of both x and f c these are quite different. Past 
(SPS) and present (RHIC) experiments suggest that the value of x may not be very 
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Figure 3: Comoving survival probability in the scenario II as a function of impact 
parameter. 

high even at LHC energies. To see the effects of x and f c on the dilepton yield we 
show in fig. (4) the invariant mass distribution for various combinations of these two 
parameters in the invariant mass range 3 GeV< M < 3.2 GeV. We notice that for 
the combination in which both x and f c are large, the yield is less, indicating more 
suppression as the beam energy increases (since in that case x will be different from 
zero). 

Now we compare the hard contribution with the thermal one in fig. (5). For 
the thermal part we use simple Bjorken cooling law. We assume a hadronic matter 
initial state with initial temperature of the order of 300 MeV at proper time n = 1 
fm. We realise that at T = 300 MeV the hadronic matter the hadronic matter may 
dissolve to QGP, however, we take such a high temperature to show that the hard 
contribution is larger than the thermal even for such a large initial temperature. 
For the hard production we choose a centrality cut of 10 %. For nuclear shadowing 
EKS [9] parametrization has been used together with CTEQ PDF [6]. We have 
included a K factor ~ 2.5 to account for higher order processes. 

As the beam energy increases the shadowing effect in hard scattering processes 
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Figure 4: Dilepton yield in Pb — Pb collisions at y/s = 6000 GeV, for various 
combinations of (x, f c ) in the second case of comoving absorption. 
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Figure 5: Dilepton yield in Pb — Pb collisions at \/s = 6000 GeV, x = 0.2 and 
f c = 0.5 in different comoving scenarios . Thermal contributions are also included. 
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Figure 6: Dilepton yield in Pb — Pb collisions at y/s = 6000 GeV, with (WS) and 
without (WOS) shadowing for two values of absorption cross-sections. 



becomes important. We have checked that at RHIC energies the lepton pair yields 
with and without shadowing are not very different. However, there is a substantial 
difference in the yields at LHC energies (see fig. (6)). 

The effects of different comoving scenarios on the dilepton yield is clearly seen 
in fig. (7) in the invariant mass window 3 GeV< M <3.2 GeV. We see that the 
second scenario of the comoving absorption gives the maximum suppression for the 
given set of parameters. Even, the thermal contribution is larger than the hard 
contribution. The yield of lepton pairs (hard part) from the experiment at LHC 
energies might lie between these yields. 



IV. Summary and discussions 

In this work we have calculated the lepton pairs yield from thermal and hard J /if) 
decays. The effects of the various comover scenarios on the survival probability 
of the heavy quarkonia, nuclear shadowing on the production cross sections of the 
heavy quarks and the centrality of the collisions have been considered. The effects 
of shadowing and comover suppressions on the dilepton spectra resulting from the 
decays of J /if) is found to be large. The dilepton yield from the thermal J /if) is 
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Figure 7: Same as fig. (5) with 3 GeV< M < 3.2 GeV. 

found to be smaller than the contributions from hard J/ip. 

The fraction of the comover is treated as a parameter here and the lepton pair 
yields from the heavy quarkonia have been evaluated for various values of this pa- 
rameter. The increase of J/ip yield from the decays of higher charmonium states 
have been neglected here, however the yield with the inclusion of such processes 
may be realized within the parameters (x, the fraction of hard component and f c , 
the fraction of the comoving hadrons) range considered here. 
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